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A Green-function theory for the dynamic spin susceptibility in the square-lattice spin-1/2 an¬ 
tiferromagnetic compass-Heisenberg model employing a generalized mean-held approximation is 
presented. The theory describes magnetic long-range order (LRO) and short-range order (SRO) at 
arbitrary temperatures. The magnetization, Neel temperature Tjv, specihc heat, and uniform static 
spin susceptibility x S'f® calculated self-consistently. As the main result, we obtain LRO at hnite 
temperatures in two dimensions, where the dependence of Tn on the compass-model interaction is 
studied. We hnd that Tjv is close to the experimental value for Ba 2 lr 04 . The effects of SRO are 
discussed in relation to the temperature dependence of X- 
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I. INTRODUCTION 

The study of the quantum compass model (CM) for 
strongly correlated transition-metal compounds with or¬ 
bital degrees of freedom and strong spin-orbit coupling 
is an active field of research (for a review, see Ref. i)- 
In particular, quantum and thermodynamic phase tran¬ 
sitions in the two-dimensional (2D) CM were studied in 
Refs. [3,i, where the directional-ordering transition of 
the 2D Ising universality class was found for the symmet¬ 
ric CM. Depending on the method of numerical compu¬ 
tation, the temperature of the phase transition for the 
quantum CM was found in the range Tc = 0.055 — 0.058 
(in the units of the CM exchange interaction) In 
Refs. @, 0 the compass-Heisenberg (CH) model was in¬ 
troduced extending the CM by the 2D antiferromagnetic 
(AF) Heisenberg model, and the ground state and ex¬ 
cited states of the CH model were analyzed. For 5d 
transition-metal compounds with a strong spin-orbit cou¬ 
pling, such as Sr 2 lr 04 and Ba 2 lr 04 with the Neel tem¬ 
peratures Tjv = 230K and Tjv = 240K, respectively, (see, 
e.g.. Ref. 0), an effective AF Heisenberg model for pseu¬ 
dospins 5 = 1/2 with the CM anisotropy was derived 
in Ref. Q. Recently we have calculated the spin-wave 
excitation spectrum and Tjv for a layered CH model by 
means of the random phase approximation (RPA) 0 (see 
also Refs. i). An important issue is the description of 
magnetic long-range order (LRO) and short-range order 
(SRO) and of the thermodynamics at arbitrary temper¬ 
atures by a theory going beyond RPA. 

In this paper we employ a generalized mean-field 
approximation (GMFA) to the 2D CH model that is 
based on the equation-of-motion method for Green func¬ 
tions. In the framework of a more general theory for 
the dynamic spin susceptibility including the self en¬ 
ergy [ifl [ill : neglect of the self-energy corresponds 

to the GMFA. This approximation, for spin-rotation in¬ 
variant (SRI) systems also named SRI Green-function 
method (RGM), has been successfully applied to several 
quantum spin systems (see, e.g.. Refs. [ia - l2q ). 


We start from the spin-1/2 CH model on the square 
lattice, 

( 1 ) 

where v = x,y,z. The nearest-neighbor (NN) ex¬ 
change interaction parameters are JL = Jij + FT, Jfj = 

Jij + Jfj = Jij, Jij = J{5rj,ri±a„ + l^rj.riiba^), 

FT = and Ff^- = FyJr^.riia,,• We assume 

J > 0 and Fa; > Fj, > 0. The symmetric formulation of 
the model © allows us to get expressions for quantities 
with indexes v = y, z from those indicated by the index 
^ = a; by cyclic permutation. 

II. THEORY OF SPIN SUSCEPTIBILITY 

To evaluate the thermodynamic quantities in the CH 
model, we calculate the dynamic spin susceptibility 
X^(a;) = -((5^|5^))a; (((...|...))ca denotes the retarded 
two-time commutator Green function 113). Using the 
equations of motion up to the second step, we ob¬ 
tain -b ((-5^|5^));^, where = 

([z5^|5^]), iS’f^ = and -5^ = [[S^^,H],H]. For 

the model o the moment is given by the exact ex¬ 
pression 

= ^[cos(qR,)(4C'o^, + JS,Cl) - JlCl - 

i 

( 2 ) 

Here, with R = me^ + ney 

denote the spin correlation functions. 

The second derivatives — 5q are approximated in the 
spirit of the scheme employed in Refs. [TH - E^ . That 
means, taking the site representation, in —Sf we decou¬ 
ple the products of three spin operators on different lat¬ 
tice sites along NN sequences {i,j,k) as 

SfSySl=af{Spl)Sf, ( 3 ) 

s^sfsy = a^,{s!sy)sj, ( 4 ) 
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where the vertex renormalization parameters af and 
0.2 are attached to NN and further-distant correlation 
functions, respectively. After some algebra, we obtain 
-5^^ = and 


with the squared spin-excitation energy 

Kf = E[M?(q) + (i-^b)^y(q)], 


( 5 ) 


( 6 ) 


where 


<(q) = - 2^KJoiCOs(qRJ], 


( 7 ) 


= oijy,jyfy^ +0^2 

+< cos(qR,,)(4 Jo"Co".) 

- cos(qRJ« Jo", Jo« + < Jo", Jo^.Co^, 

+a^ Jo^^Jo^.C^- + Jo^, Jo",C^). (8) 

The appearance of AF LRO at T < T/v is reflected 
in our theory by the divergence of Xq = Xq(‘^ = 0) 
corresponding to the closure of the spectrum gap at the 
AF ordering vector Q = (tTjTt), Wq(T < T/v) = 0. In 
the LRO phase the correlation functions = C'^ „ are 
written as [l^, [H, [l^, [12,112, 


proceed as follows. As in our recent study of the model 
o by means of the RPA and linear spin-wave theory 
(LSWT) Q, we consider an anisotropic CM interaction. 
Lx > Lj, > 0, so that the LRO phase is an easy-axis AF 
with the magnetization along the x axis. Accordingly, we 
put = 0, where we can also consider the limit¬ 

ing case Fx = Fy. Then , at T = 0, for determining the 
seven quantities 2 C", besides the three sum rules 

and the three LRO conditions, we need an additional 
condition. To this end, we adjust the ground-state mag¬ 
netization m"(0) to the expression obtained in LSWT, 
to"(0) = mls^j.{Q) = (1/2) - {l/N) Eq(5'q 5'+), where 
the correlation function is given by Eq. (17) of Ref. @ 
with the sublattice magnetization a substituted by the 
spin S = 1/2. In the LRO phase, 0 < T < T/v, 
we have found that the ansatz fa{T) = r„(0), where 
7a = with T^iJ) = {a'i{T) - \)/[o'i{T) - 1), 

is a reasonable approximation, as will be demonstrated 
in the Appendix for the Heisenberg limit. At T > T/v, for 
calculating a'{ 2 we use the sum rules and the condition 
r-(((r) = r)((T/v). 

Because, as an input, we take the ground-state magne¬ 
tization in LSWT that describes the LRO quite well for 
small enough CM interaction Fx^^y as compared to the 
exchange interaction J, we present results in the region 
Lx.y/J < 1- 


III. RESULTS 


^ E ( 9 ) 

qJQ 

with Cq calculated from the Green function ([5]) by the 
spectral theorem, 

TH^ 

C'q = (*5;i5q) = ^[l + 2n«)], (10) 

where n{uj) = (e“/^ — 1)“^. The condensation part C' 
arising from Wq = 0 determines the staggered magneti¬ 
zation that is defined by 

= = (11) 

R 

In the paramagnetic phase, we have Wq > 0 and = 0. 
The NN correlation functions are related to the internal 
energy u per site, u = (l/2Af) 1™™ which 

the specific heat Cy = du/dT may be obtained. 

To calculate the thermodynamic properties, the corre¬ 
lation functions C^, the vertex parameters a/ 2 appear¬ 
ing in the spectrum Wq, and the condensation term C'^ 
in the LRO phase have to be determined. Besides Eqs. 
(ED and (I10|) for calculating the correlators, we have the 
sum rules = 1/4 and the LRO conditions uJq = 0; 

that is, we have more parameters than equations. To 
obtain a closed system of self-consistency equations, we 


As described in Sec. II, the thermodynamic quantities 
are calculated from the numerical solution of a coupled 
system of nonlinear algebraic self-consistency equations. 

First we consider our results on the magnetic LRO. In 
the Heisenberg limit we have no LRO at finite temper¬ 
ature, in accord with the Mermin-Wagner theorem. In 
Fig. [T]the Neel temperature T/v as function of Fy for two 
ratios of Fx/Fy is plotted. Our results for T/v remarkably 
deviate from those found in RPA Q. For Fx/Fy = 1.5, 
T/v exhibits qualitatively the same dependence on Fy, but 
is appreciably reduced as compared to RPA. For the sym¬ 
metric CM interaction Fx = Fy = F, our theory yields 
a finite Neel temperature that increases with F, whereas 
in RPA, T/v = 0 was found. We ascribe these differences 
to a better description of strong spin fluctuations by our 
Green-function theory going one step beyond RPA. Note 
that in the RPA Q the fluctuations of the x component 
of spin are not taken into account. Generally speaking, 
from our results we conclude that the CM interaction 
added to the Heisenberg model favors magnetic LRO at 
finite T in two dimensions. 

Comparing our findings to experiment, we consider the 
compound Ba 2 lr 04 , where quantum chemistry calcula¬ 
tions [111 yield the parameters J = 65meV, F = 3.4meV 
= 0.0523J, and a very small interplane coupling J^ Ri 
5 X 10“^ J that is not taken into account here. In Fig. [5] 
the magnetization m = m" as function of T is plotted. 
At T = 0, m is enhanced by F as compared to the value in 
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FIG. 1: Neel temperature Tjv for Fa; = Fy (solid) and Fa; = 
l.SFy (dashed) compared with the RPA result (dotted). 



T/J 

FIG. 2: Maenetization m and specific heat Cv for Fa, = F,, = 
0.0523J. 


the Heisenberg limit (tolsvft(O) = 0.303). The second- 
order phase transition occurs at T)v = 0.392J = 295K 
which agrees rather well with the experimental value 
rj^^p _ 240K. At Tjv the specific heat (inset to Fig. ^ 
reveals a cusplike singularity. Such a structure at the 
transition temperature was also found in layered mag¬ 
nets treated by RGM in Ref. [l^. From the analysis 
given there we suggest that also here the height of the 
cusp may be underestimated. It is desirable to compare 
our result for Cv to experimental data which, however, 
are not yet available. 

Finally, we discuss the uniform static spin suscepti¬ 
bility x(r) = Xq=o(w = 0) depicted in Fig. [3] for the 
symmetric CM interaction, which reflects the behavior 
of magnetic SRO. Considering the Heisenberg limit, the 
increase of y with temperature is caused by the decrease 
of AF SRO, i.e., of the spin stiffness against the orien¬ 
tation along a homogeneous magnetic field. Since the 
SRO is less pronounced at higher T, y reveals a max¬ 
imum near the exchange energy and a crossover to the 
high-temperature Curie-Weiss behavior. In the param¬ 
agnetic phase for F > 0, the susceptibilities x'^ low- 



T/J 

FIG. 3: Uniform static susceptibility x for Fj, = Fy = 
0.0523J, (solid), (dashed), (dotted), and for the 
Heisenberg limit (dot-dashed). 


ered as compared with the Heisenberg limit due to the F- 
induced enhancement of AF SRO. Moreover, the suscep¬ 
tibility y“ = y^ is reduced as compared with y^, because 
the SRO of the x- and y-spin components is more pro¬ 
nounced due to the AF CM interaction along the x- and 
y-axes. In the LRO phase, the susceptibility strongly 
decreases with decreasing T, because the increasing easy- 
x-axis LRO also enhances the SRO. Considering the sus¬ 
ceptibility y*^, with decreasing T below Tjv the SRO of 
the y-spin components is reduced due to the increasing 
alignment of the spins in x-direction, so that x^ increases 
with decreasing T. Thus, y^ exhibits a minimum at Tjy. 
The susceptibility y^ in the LRO phase turns out to be 
nearly temperature independent. 


IV. SUMMARY 

In this paper we have calculated thermodynamic 
properties of the 2D antiferromagnetic CH model 
within a generalized mean-field approach for arbitrary 
temperatures. Our main focus was the detailed investi¬ 
gation of magnetic LRO. We have found that the Neel 
temperature Tjv is finite even in the case of a symmetric 
CM interaction. We have explained the temperature 
dependence of the uniform static susceptibility in terms 
of magnetic SRO. Our investigation forms the basis 
for forthcoming extended studies of the CH model. 
Especially, the generalization of the theory holding 
for arbitrary CH model parameters and allowing for 
the description of magnetic and directional-ordering 
transitions is of particular interest. 

The authors would like to thank W. Janke and J. 
Richter for valuable discussions. The financial support 
by the Heisenberg-Landau program of JINR is acknowl¬ 
edged. 
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TABLE I: Heisenberg limit at T = 0. 



V = X 

v = y,z 

average 

SRI 

Cio 

-0.1542 

-0.0989 

-0.1173 

-0.1173 

C 20 

0.118 

0.0412 

0.0668 

0.0667 

Cii 

0.125 

0.0522 

0.0765 

0.0763 

Jx 

0.0327 

0.0636 

0.0533 

0.05292 

n/j 

3.477 

2.821 

3.04 

2.978 


Appendix A: Heisenberg limit 

In the limit = 0, the LRO in the ground state of 
the Heisenberg AF may be described within the RGM 
by the SRI forms of Eqs. Q and (fTTl) . i.e., by = 
= C'n and as it was done, e.g., in 

Refs. (T^fisl According to the consideration 

of the CH model with the easy-axis magnetization 
let us outline an alternative possibility to describe the 
ground-state LRO in the AF Heisenberg model. Here, 
we may also break the rotational symmetry by putting 
= 0, which is analogous to the introduction of 
a symmetry-breaking sublattice magnetization (see, e.g.. 


our RPA approach Q). Performing the calculations at 
T = 0 as described in Sec. H, we obtain the correlation 
functions and the uniform static spin susceptibility 
X'^ listed in Table HI The spin-excitation spectrum at T = 
0 is found to be 

< = (Al) 

where 7 q = ^(cos^a, -I- cosqy) and = 

— 16J‘^af{CfQ + Cfg). The amplitudes ff' are given in 
Table ffl^Note that Eq. (lAip leads to the SRI result 
(Ref. ^5[), if we put = ai and C^q + C^q = 2Ciq, 
and has the same shape as in LSWT and RPA. 

To relate the non-SRI description of LRO to the 
SRI formulation, we calculate arithmetic averages, e.g., 
Cmn = Can be seen from Table U we 

obtain a very good agreement of the averaged ground- 
state properties with those resulting from the SRI the¬ 
ory. Moreover, the geometrical average fQ(0) used in 
our ansatz for the LRO phase (see Sec. H) given in the 
Heisenberg limit by fQ(0) = 1.1913 nearly agrees with 
the ratio rft(O) = 1.2109 obtained in the SRI approach 
(see also Ref. [l^)- This gives some justification for for¬ 
mulating our ansatz in terms of the geometrical average. 
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